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Abstract. In this paper I examine the phase space dynamics in the framework 
of Non- Project able Hofava-Lifshitz bouncing cosmologies. By considering a closed 
Friedmann-Lemaitre-Robertson- Walker (FLRW) geometry, the first integral contains 
a correction term that leads to nonsingular metastable bounces in the early evolution 
of the universe. The matter content of the model is a massive conformally coupled 
scalar field, dust and radiation. A nonvanishing cosmological constant connected to 
a de Sitter attractor in the phase space is also assumed. In narrow windows of the 
parameter space, labeled by an integer n > 2, nonlinear resonance phenomena may 
destroy the KAM tori that trap the scalar field, leading to an exit to the de Sitter 
attractor. As a consequence nonlinear resonance imposes constraints on the parameters 
and in the initial configurations of the models so that an accelerated expansion may 
be realized. 
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1. Introduction 

Although General Relativity is the most successful theory that currently describes 
gravitation, it presents some intrinsic crucial problems when one tries to construct a 
cosmological model in accordance with observational data. In cosmology, the ACDM 
model gives us important predictions concerning the evolution of the universe and about 
its current state [1]. However, let us assume that the initial conditions of our Universe 
were fixed when the early universe emerged from the semi-Planckian regime and started 
its classical expansion. Evolving back such initial conditions using the Einstein field 
equations, we see that our universe is driven towards an initial singularity where the 
classical regime is no longer valid [2j. 

Notwithstanding the cosmic censorship conjecture [3], there is no doubt that 
General Relativity must be properly corrected or even replaced by a completely new 
theory, let us say a quantum theory of gravity. This demand is in order to solve the 
issue of the presence of the initial singularity predicted by classical General Relativity 
in the beginning of the universe. 

One of the most important characteristics of our Universe supported by 
observational data is its large scale of homogeneity and isotropy. However, when 
we consider a homogeneous and isotropic model filled with baryonic matter, we find 
several difficulties by taking into account the primordial state of our Universe. Among 
such difficulties, we can mention the horizon and flatness problems [1]. Although the 
Inflationary Paradigm[l] allows one to solve problems like these, inflationary cosmology 
does not solve the problem of the initial singularity. 

On the other hand, since 1998 [5] observational data have been giving support to the 
highly unexpected assumption that our Universe is currently in a state of accelerated 
expansion. In order to explain this state of late-time acceleration, cosmologists have 
been considering the existence of some field - known as dark energy - that violates 
the strong energy condition. Although it poses a problem to quantum field theory 
on how to accommodate its observed value with vacuum energy calculations [6], the 
cosmological constant seems to be the simplest and most appeahng candidate for dark 
energy. Therefore, nonsingular models which provide late-time acceleration should be 
strongly considered. 

During the last decades, bouncing models [3 [8] have been considered in order to 
solve the problem of initial singularity predicted by General Relativity. Such models 
(as [9]) might provide attractive alternatives to the infiationary paradigm once they can 
solve the horizon and fiatness problems, and justify the power spectrum of primordial 
cosmological perturbations inferred by observations. 

In 2009, P. Horava proposed a modified gravity theory by considering a Lifshitz- 
type anisotropic scaling between space and time at high energies [lO] . In this context , it 
has been shown [11], [1^ that higher spatial curvature terms can lead to regular bounce 
solutions in the early universe. Since its proposal, several versions of Hofava-Lifshitz 
gravity have emerged. 
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In the case of a 4-dimensional (1 + 3) spacetime, the basic assumption which is 
required by all the versions of Horava-Lifshitz theories is that a preferred foliation of 
spacetime is a priori imposed. Therefore it is natural to work with the Arnowitt-Deser- 
Misner (ADM) decomposition of spacetime 

rfs^ = N^dt'^ - QijiN'dt + dx'){N^dt + dx^), (1) 

where N = N{t, x*) is the lapse function, A^* = N^{t, x*) is the shift and Qij = gij{t, x*) 
is the spatial geometry. In this case the final action of the theory will not be invariant 
under diffeomorphisms as in General Relativity. Nevertheless, an invariant foliation 
preserving diffeomorphisms can be assumed. This is achieved if the action is invariant 
under the symmetry of time reparametrization together with time-dependent spatial 
diffeomorphisms. That is: 

t^t{t), x' ^x'{t,x'). (2) 

It turns out that the only covariant object under spatial diffeomorphisms that contains 
one time derivative of the spatial metric is the extrinsic curvature Kij 



2N 



(3) 



L dt 

where Vj is the covariant derivative built with the spatial metric gij. Thus, to construct 
the general theory which is of second order in time derivatives, one needs to consider 
the quadratic terms KijK^^ and - where K is the trace of Kij - in the extrinsic 
curvature. By taking these terms into account we obtain the following general action 

5 oc y Ny/^[KijK'^ - - U{g^j, N)]dx^dt (4) 

where g is the determinant of the spatial metric and A is a constant which corresponds 
to a dimensionless running coupling. As in General Relativity the term KijK"^^ — 
is invariant under four-dimensional diffeomorphisms, we expect to recover the classical 
regime for A — )■ 1. That is why it is a consensus that A must be a parameter sufficiently 
close to 1. In general, U{gij, N) can depend on the spatial metric and the lapse function 
because of the symmetry of the theory. It is obvious that there are several invariant 
terms that one could include in U. Particular choices resulted in different versions of 
Hofava-Lifshitz gravity. 

Motivated by condensed matter systems, P. Hofava proposed a symmetry on U 
that substantially reduces the number of invariants |10j. In this case, U depends on a 
superpotential W given by the Chern-Simons term, the curvature scalar and a term 
which mimics the cosmological constant. It has been shown [13] that this original 
assumption has to be broken if one intends to build a theory in agreement to current 
observations. 

The simplification A^ = N{t) was also originally proposed by Hofava [lOj. This 
condition defines a version of Hofava-Lifshitz gravity called Projectable. As dN/ dx^ = 0, 
the Projectable version also reduces the number of invariants that one can include in U. 
The linearization of this version assuming a Minkowski background provides an extra 
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scalar degree of freedom which is classically unstable in the IR when A>lorA<l/3, 
and is a ghost when 1/3 < A < 1 |13]. Although some physicists argue that higher 
order derivatives can cut off these instabilities, it has been shown [T3l flSl fT6l fT7] that a 
perturbative analysis is not consistent when A — )• 1 and the scalar mode gets strongly 
coupled. That is because the strongly coupled scale is unacceptably low. In this case, 
higher order operators would modify the graviton dynamics at very low energies, being 
in conflict with current observations. 

Besides pure curvature invariants of Qij, one may also include invariant contractions 
of d{\n N)/dx'^ in U. This assumption defines the so-called Non-Projectable version of 
Hofava-Lifshitz gravity. Connected to the lowest order invariant di{ln N)d^ [In N) , there 
is a parameter a which defines a "safe" domain of the theory [18] . In fact, in this 
case there is also an extra scalar degree of freedom when one linearizes the theory in a 
Minkowski background. However, when < a < 2 and A > 1 this mode is not a ghost 
nor classically unstable (as long as detailed balance is not imposed) . Although the Non- 
Projectable version also exhibits a strong couphng[13l[T8l[19], it has been argued that its 
scale is too high to be phenomenologically accessible from gravitational experiments |14j. 

In this paper I adhere to the so-called Non-Projectable Hofava-Lifshitz gravity in 
which I consider a nonsingular FLRW cosmological model [IS]. The matter content 
is given by dust, radiation and a conformally coupled scalar field. I also assume a 
nonvanishing cosmological constant in order to obtain a de Sitter atractor in the phase 
space. In this context I show how an alternative exit to late-time acceleration (connected 
to the de Sitter attractor) may be realized. 

In the next section I present a nonsingular homogeneous and isotropic cosmological 
model - sourced with perfect fluids, a cosmological constant, and a conformally coupled 
scalar field - which arises from Non-Projectable Hofava-Lifshitz gravity. In section 3 
I analyze the structure of the phase space. In section 4 I restrict myself to the case 
of dust and radiation in order to construct a simple model. In section 5 I show how 
nonlinear resonance can provide an exit to the de Sitter attractor. In section 6 I exhibit 
the pattern of the resonance windows and show in which regions in the parametric space 
late-time acceleration may be realized. Final remarks are given in the Conclusions. 

2. The Model 

Let us consider a model in which the matter content is given by a nonminimally coupled 
massive scalar field and noninteracting perfect fluids with equation of state pi = UiPi 
(i=l,.., N). In this context, the 4-D covariant Lagrangian Cm for the matter content can 
be written as 



TV 




(5) 



where £j are the Lagrangians for the noninteracting perfect fluids and 




(6) 
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with R being the 4-D Ricci scalar. That is, is the Lagrangian density of the massive 
scalar field plus perfect fiuids whose dynamics interact only with the metric Qa/s- We 
further assume that the scalar field is nonminimally coupled with gap, with coupling 
parameter ^. 

The fundamental symmetry assumed in Hofava-Lifshitz gravity (invariant under 
diffeomorphisms that preserve the foliation) provides enough gauge freedom to choose 

N = l, N' = 0, g,^ = -a{t)\,. (7) 

This puts the geometry ([T]) into the FLRW form 



ds^ = -dr + a{ty 



(8) 



1 — kr^ 

where k is the spatial curvature, a{t) is the scale factor, t is the cosmological time and 
(r, 6, x) are comoving coordinates. It's straightforward to show that the energy density 
connected to Cm is given by 



H + 



ki 



(9) 



with the d = da/dt and H = d/a. 

By considering the local Hamiltonian constraint in the Non-Projectable version of 
Hofava-Lifshitz gravity [18], we obtain the following first integral 



2 ja^k^ a2k^ k A ^ 47rG' 
y ^ 3A - 1 1 6a4 ^ 2a2 ^ 2 ~ 



0. 



(10) 



where a2 and are constants coupling coefficients. From ( jlOl) we notice that the 
correction term proportional to a2 behaves just like a radiation component. 

Let T^^ be the energy momentum tensor connected to the matter content connected 
to Cm- As in this case the conservation equations VqT"'^ = still apply, we obtain 

k + SH{p,+p,) = 0-,p, = -^ (11) 

for the noninteracting perfect fiuids, where Ei are constants of motion. On the other 
hand, for the nonminimal coupled scalar field we obtain the following equation of motion 



1 



J 



La 



0. 



(12) 



If < 1 for all i, the conditions for the bounce are given by 

A > 1/3 , ask^ > 0. (13) 

From now on I will restrict myself to the case of closed geometries, that is = 1. It 
is worth to mention that this model is not the only possibility to generate a bounce 
in non-relativistic theories. In fact, it has been shown - for the case closed of FLRW 
geometries - that the Universe can undergo through a bounce as long as the terms which 
violate relativity lead to a dark radiation component with negative energy density |20j . 
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In order to simplify the numerical analysis I will also fix A = 1 and 0:3 > 0. In the 
framework of Hofava-Lifshitz gravity, A must be close to unity in order to assure that 
severe Lorentz violation do not to occur. Thus, the assumption of A = 1 enables the 
following model to be a fair approximation derived from of a healthy non-projectable 
Hofava-Lifshitz theory. On the other hand, as > is the necessary and sufficient 
condition for the bounce. 

By choosing the so-called conformal coupling = 1/6), equations (fTOl) and ([T2|) 
may be rewritten using the conformal time rj = J ^dt as 

3a'' + V{a) - i[v9'' + (1 + m^a^)^^] = 87rG{Erad - ^4), (14) 

and 

ip" + {l + m'^a^)ip = 0, (15) 

where Erad corresponds to a constant connected to Urad = I, the primes denote 
derivatives with respect to conformal time, (p = a^v^SvrG, 0:4 = 3a2/8iTG and 

Via) = 5 + 3a^ - Aa^ - 8nG (16) 

ij^rad 

It is worth noting that for m = the system of equations f|T4l) and f|T5|) is separable 
and, therefore, integrable. That is, in this case, equation ( |T5l) has a first integral 

is a constant of motion and, from (1141) . we obtain 




I 'l87rG{Erad-a,)+S^^-V{a)'^''- ^^^^ 
3. The Structure of the Phase Space 

By considering equations f ll4p and (1151) . it can be defined the following dynamical system: 
^' = P^, (18) 
a' = ^ , (19) 
p'^ = -(A; + mV)¥., (20) 



Pa 



,2„2 

6a + 4Aa^ + ^ + 8vrG [ -^(1 - Swi) + m^aip'^] . (21) 



Eqs. f|T8|) and (19) are mere redefinitions. On the other hand, {a,Pa) can be shown to be 
canonically conjugated by considering the first integral (HM as a Hamiltonian constraint. 
Now we focus on three basic structures that organize the dynamics in the phase space 
of the above dynamical system. 

3.1. Invariant Plane 

Let us consider the arbitrary dynamical system with n degrees of freedom {ipi, ..,ipn) G 
M", whose differential equations are given by 

^ = F,(^i,..,^„) , (22) 
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where i = 1, ..,n. Let V'ilO be its solution for given initial conditions "1^,(^0) ^ -M C M". 
Thus, Ai is defined as an invariant manifold if the solution tpi{C,) G Ai, for each 
^.(^o) e M. 

If we fix the initial conditions p^o = = ipo we see from equations (fT8|) -(21) that 
the dynamics is integrable and does not evolve in the (p and ip' directions. That is, 
orbits with these initial conditions remain contained in the plane {a,pa) during all the 
evolution of the system. Therefore the invariant plane is defined by 

ip = = p^. (23) 

It is worth noting that the dynamics in this plane is analogous to that of the 
dynamics in the separable case m = 0. In fact, in both cases the dynamics is separable 
and integrable, and its description is given by similar orbits which differ by a constant 
in the {a,Pa) sector. As we shall see, in order to furnish an exit to a de Sitter attractor 
due to parametric resonance I will always assume initial conditions sufficiently close to 
the invariant plane. 



3.2. Critical Points 



In the phase space {^pl, .., ^pn) of an arbitrary dynamical system like fl22|) . a critical point 
ipi"^ is defined as a solution of the equations Fi{ipl^, .., ip^) = 0. That is, it is a stationary 
solution of ( I22I) . If one takes the initial condition ipii^^o) = V^f^, then ipiiO = i'i'^ 

e 

The structure of the phase space of (IT8l)-(21) allows the presence of critical points 



[if = 0,a 



V'iarr) = — 



0,P« 
dV 
da 



0), where a^r satisfies the relation 
= 0. 



(24) 



It's easy to see that, by definition, the critical points are contained in the invariant 
plane. Furthermore, according to fl2^ the critical points are associated to potential 
extrema. For specific numerical values of A, Ei and Ui, we may obtain one or many 
extrema for V{a) (characterized by one or many values of Ocr). In fact, that is the case 
for a fixed value of A and suitable domains of E^. For A = 0, the dynamical system 
(fT8l)-(21) has only one critical point connected to a global minimum of the potential 
V{a). As an exit to the de Sitter attractor can not be performed in this case, I will not 
consider such configurations. 

Linearizing the dynamical system fl22l) around the critical points we obtain 

dipi 
d^ 



(25) 



where 



dFi 
dtpi 



dFi 



(26) 
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Let us then assume that the matrix dFij has fii = idi pure imaginary eigenvalues (where 
/ = 1, .., 2p, p being an integer smaller than n/2) and = Xg (where s = 2p+l, ..,n) real 
eigenvalues. The respective eigenvectors are given by x; = w / + iq ; and v^. Therefore 
we define the following local subspaces (in a small neighbourhood of the critical point) 
in the phase space 

W' = Span{vs\Xi < 0}, (27) 

W = Span{vs\Xi > 0}, (28) 

W'' = Span{wi,cii\Xi = 0}. (29) 

The superscript s, u and c denote stable, unstable and center manifolds respectively. If 
dFij has only pure imaginary eigenvalues, the critical point is called a center. If dFij 
has n — 2 pure imaginary eigenvalues and two real eigenvalues (one positive and one 
negative), the critical point is called a saddle- center. 

In the case of the dynamical system ( fT8i) -(21). we obtain 

( ^' \ ( 1 

-m^a^J 




a' 

\p'a I 



\ 











da? 



\ 


/ 




\ 




1 

6 




Pv 




a 


- Q-cr 




0/ 


\ 


Pa 


1 



(30) 



It can be shown that the above matrix has the following eigenvalues 

±1 



±i 



1 + vn?al^^^ , /i3,4 



6 



(31) 



By considering the plane {ip,p^), the corresponding eigenvectors of /ii and /i2 engender 
the local topology of around the origin {ip = 0,p^ = 0). Thus, the local topology of 
the critical points is determined by the second derivative of the potential V{a). 



When 



da'2 



> we obtain a local minimum for the potential V{a) and one 

—0,cr 

more pair of pure imaginary eigenvalues. In this case, by considering the plane (a,pa), 
the corresponding eigenvectors of /is and fi^ engender the local topology of around 
,Pa = 0). As a consequence, the topology around the critical points 
> is given by 5*^ x 5*^ (or 2-torus). We denote by Pq a critical 



the point (a 
with 



point with this property. It is called a center by definition. 



When 



9a2 



< we obtain a local maximum for the potential V{a) and 

—O'er 

two real eigenvalues (one positive and one negative). In this case, by considering the 
plane (a,pa), the corresponding eigenvectors of /is and /i4 engender the local topology of 
a saddle around the point ( = 0). As a consequence, the topology around the 

critical points with ( ^ gjt'' ) < is given hj R x S^. We denote by Pi a critical 

\ / a=acr 

point with this property. It is called a saddle-center by definition. 

The expansion of the first integral ( IT^ around the critical points reads 



H 



Y^pI + ^^"(«c«i) (a 



-87rG(E,,rf - «4) + 0{3) 



0, 



-\[pl + {'^+^^O.litW]+E,r^t 



(32) 
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where 0{3) denote terms of higher order in the expansion and Ecru = V{0'crit) is the 
energy of the respective critical point. In a small neighborhood of the critical points 
these higher order terms can be neglected and the dynamics is nearly separable in the 
sectors (a,Pa) and {f,Pip) with constants of motion given by 

1 2 1- 
1 



^{«) = TT^pI + o^"(^crit) (a - Ctcrit)^ , (33) 



p'+il + m'ai-M , (34) 



with E(^a) - -E(^) + Ecrit - 87iG{Erad - ^4) ~ and \Ecrit - SnGiErad - "4)1 sufficiently 
small. While in the sector {(f,p^) we have a rotational motion with energy around 
the critical points, in the sector {a,Pa) we have two possibilities. If V"{acr) > we have 
a rotational motion with energy E(^a) in a small neighborhood of the critical point Pq. 
Otherwise (y"{acr) < 0), we obtain a hyperbolic motion around the critical point Pi. 
The critical point Pi defines an universe analogous to that of the unstable Einstein static 
universe [2]. On the other hand, the configuration of stable Einstein static universe, 
corresponding to the critical point Pq, does not possess any classical analogue. 



3.3. Separatrices 

According to the definition above, from the saddle-center emerges a special structure 
consisting in two local subspaces. While one of them is generated by the unstable 
manifold W^, the other is generated by the stable manifold W^. Being transversal to 
each other they define the separatrices of the saddle- center critical point. 

Let us then assume that one of the initial conditions is given by a point in the 
phase space which lies on W"^ in a neighbourhood of the saddle-center critical point. 
Although this manifold is locally unstable, in general this does not mean that the final 
stage of the dynamics differs from the saddle- center critical point. In fact, sometimes 
the nonlinearities of the dynamics may induce an orbit to join the saddle-center point 
to itself. In this case, the final stage of the dynamics is the the very same saddle-center 
critical point. Orbits with such a property are called homoclinic orbits. 

From the saddle-center critical point Pi (when present) emerges a structure of 
separatrices S contained in the invariant plane. One of them tends to a de Sitter 
attractor at infinity, defining an escape of orbits to an accelerated phase regime. In 
fact, a straightforward analysis of the infinity of the phase space shows the presence of 
a pair of critical points in this region, one acting as an attractor (stable de Sitter 
configuration) and the other as a repeller (unstable de Sitter configuration). The 
scale factor approaches the de Sitter attractor as a{ri) ~ (Co — rj)^^ for rj — )• Co, or 
a{t) ~ exp(t^/A/3). 
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4. A Simple Model 

Let us now consider that the noninteracting perfect fluids are given by dust and 
radiation. In this case it can be shown that the potential 

V(a) = 3a^ - Aa^ - MGEdusta + ^ (35) 

will have two extrema (one local minimum and one local maximum) as long as the 
following conditions hold: 

21{MGEaustf > (36) 



and 



where 



6^ - 473A - 8nGEd^,t > (37) 
^ 12A ^ p ' 



256 

lOSiSnGEdust) + I2V3 J 27i8TTG Edust)^ 



A 

It can be numerically shown (cf. Fig. 1) that the increase Edust has the effect of spoiling 
the potential conflguration with two extrema. 

If m = the dynamics is integrable and thus separable. In this case the flrst integral 
( IT41) reads 

2 

^ + Via) = 8nG{Erad - a,) + £j = (38) 

and the scalar fleld behaves just like a radiation component in the dynamics of the scale 
factor. 

In Fig. 2 I exhibit the phase portrait in the invariant plane = = p^. With 
suitable values for the parameters, the two extrema of V{a) are connected to Pq and Pi. 
This model furnishes us with perpetually bouncing universes (periodic orbits in region 
I). The two separatrices Si and S2 that emerge from Pi coalesce generating a boundary 
for region I. This boundary deflnes an homoclinic orbit by definition. Orbits in region 
II and III correspond to universes with one bounce only. 

Now I focus on some structural differences between the integrable dynamics in the 
invariant plane and the integrable dynamics given by m = 0. If (p{0) and/or p^(0) do 
not vanish (in the integrable case m = 0), the phase portrait in the plane {a,Pa) is 
analogous to that of the invariant plane (cf. Fig. 2). Let Xq and Xi be the analogous 
that of Pq and Pi (connected to the potential extrema). In this case, Xq and Xi would 
be no longer critical points. Instead, they define stable and unstable periodic orbits 
respectively. Due to the system's periodicity in a neighborhood of Pq, orbits are confined 
on two dimensional surfaces in the phase space which topologically coincide with 2-tori. 
Therefore, the integrable dynamics (m = 0) is not constrained in the invariant plane 
but, in a neighborhood of Pq, evolves on 2-tori which are the direct product of closed 
curves (analogous to that of in region I) with periodic orbits in the sector {ip,p^). 
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Figure 1. The potential V{a) for several values 8 = E^ust- For higher values of 5 the 
minimum of the potential is no longer present. For 5 ~ 2.3 the two extrema of the 
potential (connected to the stable and unstable Einstein static universes) vanish. Here 
we fixed A = 1.5, as = 10~^ and SttG = 1. 




Figure 2. The phase portrait of invariant plane dynamics with the critical points 
Po (center) and Pi (saddle-center) corresponding to stable and unstable Einstein 
universes. The periodic orbits of region I describe perpetually bouncing universes. 
Orbits in Region II and III are solutions of one-bounce universes. A separatrix 
emerges from Pi defining an escape to the dc Sitter attractor. Here we fixed A = 1.5, 
as = 10-3 = Edust, Erad = 0.5 and SttG = 1. 
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Another topological feature to be considered - in the integrable case m = - is the 
dynamics of orbits in a neighborhood of Pi. Given the initial conditions Pa{0) = and 
flo = cl\xi, the motion of the system corresponds to an unstable periodic orbit which 
generates a circle in the sector {ip,p^). Let T denotes one of these orbits. Therefore, 
the direct product of periodic orbits in the sector (v2,p^) with the stable and unstable 
manifolds (analogous to that of 5*1 and S2) generates semi-cylinders which coalesce in 
the very same orbit F. When we consider a nonlocal description of the topology, the 
nonlinearities induce such cylinders to close onto themselves. As the projection of such 
cylinders in the sector {a,Pa) defines an homoclinic orbit, we call these structures as 
homoclinic cylinders. On the other hand, the direct product of periodic orbits in the 
sector {(p,p,p) with the analogous to that of separatrices 5*3 and 5*4 generates semi-infinite 
cylinders which coalesce in the orbit F. A summary of structural differences between the 
dynamics in the invariant plane and the integrable case m = is given in the table below. 



Dynamics in the Invariant Plane: 


Integrable Dynamics m = 0: 


No motion in the sector {ip,p^) 


Separable motion in the sectors {a,Pa) and {(p,Ptp) 


Critical points Pq e Pi 


Stable and unstable periodic orbits 


Periodic orbits (Region I) 


Integrable Tori 


Separatrices between regions I e II 


Homoclinic cylinders 


Separatrices between regions II e III 


Semi-infinity integrable cylinders 



5. Nonlinear Resonance of KAM Tori 

According to Liouville- Arnold [21] theorem, if the motion of a hamiltonian system with 
n degrees of freedom is integrable and bounded, then the orbits of such a system are 
confined on n-dimensional hyper-surfaces in the phase space which topologically coincide 
with n-tori. That is exactly what occurs with our system when we consider the integrable 
case m = in a neighborhood of Pq. 

In order to give a more precise analysis, let us consider the surfaces with energy 
= SnG^Erad — 0^4) + ° < Ecrit{Pi) = V{Pi). This region of the phase space is 
foliated by 2-tori x which are the topological product of periodic orbits of the 
separable sectors {^,Pip) and (a,Pa). £° and are conserved quantities for those orbits. 
The frequency Va of the periodic orbit in the sector (a,Pa) is given by the integral 

L.Mr ^ . (39) 

where /3j {i = 1...6) are the real roots of V{a) = 87iG{Erad — ^4) + Here I denote /3j 
(i = 1...3) as the three positive real roots with Pi < < (3^. On the other hand, the 
periodic orbits in the sector ((^,p^), parameterized by have the frequency i/^ = l/27r. 
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The importance of n-tori in integrable hamiltonian systems comes from the fact 
that these surfaces trap the dynamics in a finite region of the phase space. In our case, 
such tori in a neighborhood of Pq avoid an exit to the de Sitter attractor. However, a 
relevant question which arises is whether such tori "survive" when one introduces a small 
perturbation connected to the mass m of the scalar field. Assuming initial conditions 
iVo^Pipo) sufficiently close to the invariant plane, equation f|T5|) may be rewritten as 



ip" + {l + m^al{r)) ip = 0, (40) 

where ao('7) is the background solution for the scale factor of the integrable dynamics 
(ITTl) . Equation ( HOl) is a Lame equation. Defining z/^^ as the frequency in the {(p,p,fi) 
given by fHOl) . the resonance phenomena |2l] will occur when the ratio 

~ ^ (41) 
is a rational number. Expanding ao('7) in the Lame equation, one can show that 



+ + (42) 



1 / (0.9 m)2 
2^V^ ^ 2" 

However, as the system evolves the amplitude of the scalar field may grow so 
that the solution of the integrable case ao{ri) is no longer a good approximation to be 
introduced in (ITSl) . As we shall see in the next section this process may lead the dynamics 
into a more unstable behavior, with the amplification of the resonance mechanism and 
the break of the KAM tori [22] . To analytically show this behavior, let us now consider 
the following approximation of constraint ([1] 



n^E'^-Sy ^m'alivVir]) ^ SnGiErad - a,), (43) 

where (f{T]) is an approximate solution of the Lame equation. Now I introduce 
the action-angle variables pT] (Oip, JZ^, 6a, J7a)- The angle variables are defined by 
(0(p = u^rj, Qa = ^aV^) in such a way that they span the interval [0, 1] during a 
complete period of the system. Taking into account that the function ao(?7) is periodic 
with period Ta = the expansion of the non- integrable term of ( H3l) is given by 



- ^m^al{r])^\7]) = ^ (a„ cos2n7re„) cos47re^, (44) 

n 

where An are constant coefficients. The superior index in J7a and J'^ denotes that these 
are the action variables for the integrable case. The Hamilton equation for J'a, derived 
from ( H3|) . can then be integrated furnishing us in first approximation with 



2 ^ ^-^ ZTmi/.r 



cos{2nnQa — 47r9^) 
R-2/n 

cos(27™9a + AttQ^) ] ^^^^ 



^ R + 2/n 
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From (1451) we see that the dominant resonance terms are those for which i? ^ 2/n. It 
can be shown that the condition n > 2 must hold in order to obtain a real positive 
numerical value for the mass m. Therefore, 



is a good approximation in order to determine the resonances of the dynamical system 
(fT8|) -(21) in the presence of dust and radiation. When such resonances occur one can 
eventually obtain a loss of stability with the break of the KAM tori [22], allowing the 
dynamics to an exit to the de Sitter attractor. 

Let us now consider a Poincare map in the variables {(p,p^) with section pa = 0. 
As a convention I will assume that this map is unidirectional. That is, a chosen orbit 
of the system crosses the plane of section Pa = only once after a period of time T^. If 
= 0, the periodic orbits in the sector {a,pa) correspond to the point (y? = 0,p^ = 0). 
For and m = 0, the tori are characterized by closed curves around the origin 

of this map. For a small value of m this configuration is maintained around the origin 
{(f = 0,p^ = 0). In fact, according to the KAM theorem [21], if Ua and are sufficiently 
irrationals (that is, satisfy the diophantine condition) then the solutions of the perturbed 
system are quasi-periodic for a sufficiently small value of m. 

In order to compare the approximation P6l) with the exact dynamics I will perform 
a numerical analysis of the evolution of the system in the Poincare map with section 
Pa = 0. For computational simplicity I will fix A = 1.5, Erad = 0.5, as = lO^'^ = Edust 
and 8ttG = 1. In this way one can define the parametric space (0:4, m) where the 
resonances may occur. For several initial conditions around {pi^^ = 0,ipo = 0) I construct 
the Poincare map with m = 9.6 (Fig. 3) and m = 11.5 (Fig. 4). According to 
approximation ( H6l) the first map shows the resonant behavior of the system for n = 3. 
The second map shows the pattern of parametric stability in a region between the 
resonances n = 3 and n = 4. The structure of the stochastic sea [25] in Fig. 3 shows 
that initial conditions near the invariant plane can generate orbits with a long time of 
diffusion before escaping to the de Sitter attractor. 

In Fig. 5, I numerically construct the resonance chart using the exact dynamics. 
Taking the initial conditions p^o = = p^^ and ipo = 10"^, the value of oq is obtained 
by substituting the numerical values of and m in the constraint (fT4l) . For a suitable 
value of ^4, the approximate expression (H6l) is an accurate guide in order to localize 
the respective values of m in which the resonances in the parametric space (04, m) 
occur. The dashed curves (cf. Fig. 5) in the parametric space {a4,m) are constructed 
using approximation (H6l) and they allow us to localize a given domain of resonance. 
As I previously pointed out, the dominant resonances of the system are connected to 
the bifurcation of stable periodic orbits at the origin. Although approximation P6l) 
allow us identify the curves in the parametric space (04, m) where the resonances occur, 
the effect of the exact dynamics tends to stretch these domains. In fact, as one may 
numerically verify, for a fixed value of 04 there is a continuum domain of values of m 
(where the bifurcation of stable periodic orbits at the origin occurs) for each resonance. 




(46) 



Nonlinear Resonance in Hofava-Lifshitz Bouncing Cosmologies 



15 









































1 ^ i; . 


























































■■■■■ ■■r^i-ij.t;^;,; 





9 

Figure 3. Poincare map with section Pa = for several initial conditions with m = 9.6 
in the domain of parametric resonance n = 3. Here we see that the resonance of the 
exact dynamics is connected to the bifurcation of the stable periodic orbit at the origin 
{ip = 0,Pip = 0). That is, when the resonance occurs this stable periodic orbit turns 
into a unstable periodic orbit. In this case no KAM tori are present around the origin 
of the map so orbits with initial conditions in a neighborhood of the origin may escape 
to the de Sitter attractor. The presence of such bifurcation is a crucial feature in order 
to allow the dynamics to an exit to the de Sitter attractor. 



0.07D 



0.035 - 
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Figure 4. Poincare map of section Pa = for several initial conditions with m = 11.5 

in the domain of parametric stability between the resonances n = 3 and n = 4. The 
topology of the system around the origin (which defines the stable periodic orbit) 
corresponds to a 2-torus. In this case, the remaining KAM tori from the integrable 
case trap the orbits with initial conditions in a neighborhood of the origin avoiding an 
escape to the de Sitter attractor. Therefore the region of parametric stability of the 
system does not favor late-time accelerating scenarios. 
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n=2 n=4 n=5 



5 10 15 20 25 30 m 

Figure 5. Resonance chart in the plane (a4,m). The dashed hnes are solutions of 
the conditions of resonance (j46|). The hatched areas correspond to resonance domains 
of the exact dynamics. These domains correspond to regions where the bifurcation 
of periodic orbits at the origin (in the Poincare map of section pa — 0) occurs. The 
white remaining region corresponds to the domain of parametric stability where the 
dynamics is trapped by the KAM tori. For computational simplicity I am still fixing 
the parameters A = 1.5, Erad — 0.5, — 10^-^ = Edust and SttG — 1. 

The windows of exact resonance are shown as hatched regions in Fig. 5. 
6. The Resonance Pattern 

The resonance regions in the parametric space (a4, m) possess a substructure which I 
now examine. In order to simphfy this analysis I will restrict myself to the resonance 
domain = 3 (cf. Fig. 5) with 04 = 0. For these fixed values the resonance occurs in 
the interval Am = [8.7,9.8]. In this interval one can notice three distinct regions. 

(i) For m < 9.3 the dynamics is highly unstable and the resonances provide a rapid 
escape to the de Sitter attractor. In Fig. 6 the behavior of a and ip (with respect to the 
conformal time) is shown by taking m = 8.8. In this figure one can observe that a rapid 
escape to the de Sitter attractor occurs when ry ~ 70 so there is no enough recurrence 
in order to construct a Poincare map. 

(ii) When 9.7 < m < 9.8 the motion of orbits is resonant and chaotic, although 
stable. In Fig. 7 the behavior of an orbit in this region is exhibited for m = 9.8. This 
is the pattern close to the right edge of resonance n = 3. Due to its stability these 
orbits are not interesting from the late-time accelerating point of view. It is worth 
mentioning that this behaviour is in agreement to the analytical description regarding 
cyclic universes given in [20] . 

(iii) A region of transition occurs when 9.3 < m < 9.7. In this case, orbits go 
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Figure 6. The evolution of the scale factor and the scalar field for m = 8.8. 
This behavior characterizes a region of disruptive resonance close to the left edge 
of resonance n = 3 of Fig. 5. 
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Figure 7. The evolution of a and Lp ior m = 9.8 (corresponding to the right edge of 
resonance zone n = 3 of Fig. 5). The motion in this region is not interesting from the 
late-time accelerating point of view. 



through a long time of diffusion before escaping to the de Sitter attractor. In Fig. 8 a 
Poincare map (with pa = 0) of an orbit with m = 9.7 is shown. This map illustrates 
what happens in the above interval. This is an example of how orbits can go through a 
long time of diffusion before escaping to the de Sitter attractor. 

The above substructure is a pattern which is maintained for every value of in the 
resonance zone n — 3. Furthermore, it can be shown that this pattern is qualitatively 
equivalent for every value of n. Throughout this analysis one notices that nonsingular 
perpetually bouncing models from Hofava-Lifshitz possess a restrict domain in the 
parametric space where late-time acceleration (connected to the de Sitter attractor) may 
be realized. For typical variations of the parameters Edust and/or 03 the domains of the 
parametric space (0:4, m) - where the system is resonant - can be stretched or shrunken. 
Nevertheless the pattern in resonance windows and its substructure are maintained as 
one may numerically verify. In this sense the pattern is said to be structurally stable. 
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Figure 8. Poincare map (with section pa = 0) of an orbit with m = 9.7. This is 
the behavior in the region of transition of resonance n = 3 of Fig 5. In this case one 
can observe that the orbit go through a long time of diffusion before escaping (when 
77 ~ 10®) to the de Sitter attractor. This map exhibits a dark region connected to the 
structure of random motion of the orbit in a stochastic sea around the KAM islands. 



7. Conclusions 

In this paper I examine the effect of parametric resonance in bouncing cosmologies 
originating from Horava-Lifshitz gravity. In this context, terms arising from foliation 
preserving diffeomorphism invariance - which breaks 4-D covariance - implement 
nonsingular bounces in the early evolution of the universe. The matter content of the 
model is given by perfect fluids, namely dust and radiation. Furthermore I also assume 
a nonvanishing cosmological constant (connected to a de Sitter attractor in the phase 
space which provides late-time acceleration) and a massive conformally coupled scalar 
field. 

By considering the case of closed geometries I obtain a potential well with a local 
minimum and a local maximum (cf. Fig. 1) respectively connected to the critical points 
Pq (center) and Pi (saddle-center) in the phase space. Assuming a conformally coupled 
scalar field, the oscillatory behavior of the dynamical system around Pq might become 
metastable when the system is driven into a resonance window of the parameter space 
- labeled by an integer n > 2. In this case I determine the physical domain of the 
parameters (cf. Fig. 5) in which the breakup of KAM tori may occur, leading the 
Universe to a late-time acceleration regime. 

It is worth mentioning that, as examined in [26], a chaotic exit to accelerated 
expansion can be also realized - in the dynamical system f|T8l) -(21) - if one assumes 
initial condition sets taken in a small neighborhood of the stable separatix Si. These 
sets possesses fractal basin boundaries connected to a code recollapse/escape leading to 
a chaotic exit to an accelerated regime. 

Although the cosmological constant poses a crucial problem to quantum field theory 
on how to match its observed value with vacuum energy calculations, the cosmological 
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constant is by far the simplest explanation for the present acceleration of the Universe. 
Indeed, the ACDM standard model assumes that there exists a cosmological constant 
which becomes dynamically important when the typical scale of the Universe has the 
size of the present Hubble radius. In this sense, the model of this paper does not exhibit 
an alternative explanation for late-time acceleration. Instead, the core of this paper is 
to examine the dynamics in the phase space of the above model, showing how to provide 
an alternative exit to late-time acceleration. 
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